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Abstract 

We study hadronic charmless two-body B decays to final states involving 
pseudoscalar and vector mesons with the QCD factorization approach, includ- 
ing the contributions from the chirally enhanced power corrections and weak 
annihilations. The CP-averaged branching ratios and CP-violating asymme- 
tries are given. Most of our results are in agreement with the present measure- 
ments, but the branching ratios of some decay channels are only marginally 
consistent with the experimental observations. Considering the large uncer- 
tainties, great advances in both experiment and theory in the near future are 
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strongly expected. 
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I. INTRODUCTION 



In the standard model (SM), CP violation is described by the phase of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix. Phenomenologically, it is clear and convenient to ex- 
plore CP violation with the well-known unitarity triangle VudV*f, + VcdVch + ^d^tfe = 0- ^ 
sense, the study of B meson decays is mainly to make enough independent measurements of 
the sides and angles of this unitarity triangle. It is thus crucial to have a clear understanding 
of exclusive hadronic charmless B decays which would give some useful information about 
and/or constraints on the unitarity triangle. Recently, several experimental groups have 
reported their latest results ||I Hl^1 , and more B decay channels will be measured with great 
precision soon. With the accumulation of the experimental data, theorists are urged to gain 
deep insight into the rare hadronic B decays, and reduce the theoretical uncertainties in 
determining the CKM parameters from experimental measurements. 

Theoretically, intense investigations of hadronic charmless two-body B decays have been 
carried out in great detail with the effective Hamiltonian. Combining operator product ex- 
pansion with the renormalization group method, the effective Hamiltonian for B decays is 
generally expressed by the products of the Wilson coefficients and dimension-6 effective op- 
erators. The Wilson coefficients can be calculated reliably by perturbation theory, they have 
been evaluated to the next-to-leading order |]T^. Thus, the main task for us is to evaluate 
the hadronic matrix elements of the effective operators. Based on the naive factorization 
(NF) hypothesis |jl8|, which is supported by the argument of color transparency [|l9l, the 
hadronic matrix elements are usually parametrized into the product of the decay constants 
and the transition form factors phenomenologically. The NF approach works well for B and 
D decays, at least it can give the proper orders of magnitude of the branching ratios for B 
and D decays. However, the NF approach is a very rough approximation, and has obvious 
shortcomings: (1) In the NF framework, the renormalization scheme- and scale-dependence 
in the hadronic matrix elements are completely missing. Its predictions would be physical 
only when the hadronic matrix elements could make compensation for the renormalization 
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scheme- and scale-dependence of the Wilson coefficients. (2) There is no direct CP vio- 
lation in hadronic B decays with the NF approach because the Wilson coefficients, decay 
constants, and form factors are all real. This indicates that "nonfactorizable" contributions, 
which account for final-state rescattering and the strong interaction phase shift, are impor- 
tant. (3) As stated in it is questionable whether "nonfactorizable" effects can be simply 
absorbed into decay constants and form factors. They will have great effects on the class 
II (02 dominant) decay modes, e.g. -B°-^p°7r°, p°?7, cur], ■ ■ -. So it is imperative to go beyond 
this level of understanding. 

Recently, M. Beneke, G. Buchalla, M. Neubert, and C. T. Sachrajda suggested a QCD 
factorization (QCDF) method for hadronic B decays in the heavy quark limit, combining the 
hard-scattering approach with power counting in 1/mf, I^Tf. The hadronic matrix elements 



{MiM2\0i\B) (here Mi denotes the recoiled meson which picks up the light spectator quark 
in the B meson, M2 is the emitted meson, and Oi are effective operators) can be computed 
from ffist principles and expressed in terms of form factors and meson light-cone distribution 
amplitudes (LCDAs) only if M2 is light or an onium. In addition, when Mi is light, there 
exist hard-scattering interactions between M2 and the spectator which do not exist in the 
NF approach. The QCDF approach shows that the "nonfactorizable" contributions are 
dominated by hard gluon exchange and therefore computable systematically. Detailed proofs 



and arguments can be found in |21]; its applications to hadronic two-body B decays can be 



found in the literature 21-28 



In recent work |2^, we calculated the branching ratios and CP asymmetries for B^PP 



(where P denotes the pseudoscalar meson) within the QCDF approach. We found that 
with appropriate parameters our results are in agreement with the current experimental 
data. This encouraged us to further investigate the exclusive decays B—>-PV (where V 
denotes the vector meson) with the QCDF method. The B~^PV decays have been carefully 
studied within generalized factorization [p0| , p9| , p0| , and some decay channels have also been 
studied with the perturbative QCD approach |^T| where the form factors are believed to be 



perturbatively calculable with the assistance of the Sudakov factors. However, it is still an 



open question whether the Sudakov factor is apphcable in B meson decays . In the QCDF 



framework, it is argued that the form factors are not fully calculable according to naive 
power counting. M. Yang and Y. Yang studied the decay modes B—>-hih2 {hi = tt,K and 



/i2 = K*,p,uj) in |25| using the QCDF approach. The differences between their work and 
ours are that we consider the chirally enhanced contributions from another type of twist-3 
LCDAs 0O- of the pseudoscalar mesons, which is crucial for gauge independence, and we also 
take the weak annihilation topologies into account. H. Cheng and K. Yang pointed out that 
the annihilation contributions could be sizeable for the decays B^ipK Sometimes, 



the annihilation contributions may even be dominant, e.g., B^^K'^K*^ . In this paper, we 
will undertake a comprehensive study of the hadronic charmless decays of B^PV within 
the QCDF framework, including the effects of weak annihilation and the chirally enhanced 
contributions. 

This paper is organized as follows. Section || is devoted to the theoretical framework. 
There we calculate the "nonfactorizable" corrections to hadronic matrix elements of the 
effective operators with the QCDF approach, including the chirally enhanced power correc- 
tions and weak annihilations. The input parameters in our calculations are given in Sec. fT|. 
The numerical values and some remarks about the CP-averaged branching ratios and CP- 
violating asymmetries for B^PV decays are arranged in Sec. and Sec. 0, respectively. 
We come to our final conclusions in Sec. |V|. 

II. THEORETICAL FRAMEWORK FOR B DECAYS 
A. The effective Hamiltonian 



The effective Hamiltonian for hadronic charmless B decays is written as 17 



10 

q=u,c k=3 
-Vt [C^^Ql^ + CsgQfig] } + H.C., (2.1) 
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where Vq = VqbV*^ (for b-^d transitions) or Vq = VqbV*g (for b—^s transitions) are CKM 
factors. Cj(/i) are Wilson coefficients; they are universal and process independent, and have 
been evaluated to the next-to-leading logarithmic order (NLO) with the perturbation theory 
and renormalization group method. We list their numerical values in the naive dimensional 
regularization (NDR) scheme at three different scales in Table |. The dimension-6 local 
operators, including tree operators Ql^Ql, QCD penguin operators Qs^Qe, electroweak 
penguin operators Qj^Qio, and magnetic penguin operators Qj-yyQsg, can be expressed 
explicitly as 



= 


{Uaba)v-A{qi3Ul3)v-A, 




{Caba)v~A {qi3C/3)v-A, 


(2.2a) 


Q^2 = 


{uJ)p)v-A{.qpUa)v-A, 


Ql = 


{Cab/3)v~A{qi3Ca)v-A, 


(2.2b) 


Q3 = 


{qaba)v~Aj2(^Wv-A, 
q' 


Q4 = 


{qpK)v-AY.(q'aq'p)v-A, 

q' 


(2.2c) 


Q5 = 


{qaba)v-A '^{q(^q'i3)v+A, 

q' 


Q6 = 


{qi3ba)v-A '^{qaq'p)v+A, 
q' 


(2.2d) 


Q7 = 


3 

<?' 


Q8 = 


3 

-iqi3ba)v-Aj2^q'(^l'l3)v+A, 
q' 


(2.2e) 


Q9 = 


3 

l^{qo^o)v-A^eq\4pqp)y-A, 
<?' 


Qio = 


3 

= 2(^pba)v-Aj2^q'(^My-^^ 
q' 


(2.2f) 


Q77 = 


g 




:^m,g,a^'^(l + 75)t^;3&/3G^., 


(2.2g) 



where q' denotes all the active quarks at the scale fj, = 0{mh), i.e., q' = u, d, s, c, b. 



B. B^PV in the QCDF framework 

When the QCDF method is applied to the decays B^PV, the hadronic matrix elements 
of the local effective operators can be written as 

{PVmB) = Ff-^{0) fdx T^,{x)^y{x) + Ar^(O) fdy Tl,{y)^p{y) 

+ Cdi fdx fdy Tl\i, X, y)<^Bmvi^)^piy)- (2-3) 
Jo Jo Jo 

Here F^^^ and A^^^ denote the form factors for B^P and B^V transitions, respectively. 
^b{C)i ^v{x)i and ^p{y) are the LCDAs of valence quark Fock states for B, vector, and 
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pseudoscalar mesons, respectively. T/'^^ denote the hard-scattering kernels, which are domi- 
nated by hard gluon exchange when the power suppressed O^AqcD/^b) terms are neglected. 
So they are calculable order by order in perturbation theory. The leading terms of come 
from the tree level and correspond to the NF approximation. The order of terms of 
can be depicted by vertex-correction diagrams Fig.l(a-d) and penguin-correction diagrams 
Fig.l(e-f). T^^ describes the hard interactions between the spectator quark and the emitted 
meson M2 when the gluon virtuality is large. Its lowest order terms are 0{as) and can 
be depicted by hard spectator scattering diagrams Fig.l(g-h). One of the most interesting 
results of the QCDF approach is that, in the heavy quark limit, the strong phases arise 
naturally from the hard-scattering kernels at the order of Og- As for the nonperturbative 
part, it is either power suppressed in l/nib or separated into the form factors and LCDAs 
of mesons. 

Because the b quark mass is not asymptotically large, the power suppression might fail 
in some cases. For instance, the contributions of operator Qq to the decay amplitudes would 
formally vanish in the strict heavy quark limit. However, it is numerically very important 
in penguin-dominated B rare decays, such as the interesting channels B^ttK, etc. This 
is because Qq is always multiplied by a formally power suppressed but chirally enhanced 
factor = rnb{mi+m2) ~ ^ (-'")' '^hcrc nil and m2 are current quark masses. Therefore 
phenomenological application of QCD factorization in B rare decays requires at least a 



consistent inclusion of chirally enhanced corrections. The readers may refer to Refs. p2|j27 
for more details. 

With the above discussions on the effective Hamiltonian of B decays Eq.( |2.1|) and the 
QCDF expressions of hadronic matrix elements Eq.( p.3|) , the decay amplitudes for B-^PV 
in the heavy quark limit can be written as 

f< 10 

A{B^PV) = ^ E Y.^,aUPy\0^\B)f. (2.4) 

V ^ p=u,c i=l 

The above {PV\Oi\B) j- are the factorized hadronic matrix elements, which have the same 
definitions as those in the NF approach. The "nonfactorizable" effects are included in the 



coefficients which are process dependent. The coefficients Oj are collected in Sec. [II C| , and 
the explicit expressions for the decay amplitudes of B—>-PV can be found in Appendix. B of 
The only differences from ours are the expressions for the coefficients a^. 
According to the arguments in [^], the contributions of weak annihilation to the decay 



amplitudes are power suppressed, and they do not appear in the QCDF formula Eq.(|2.3|). 
But, as emphasized in [0,^, the contributions from weak annihilation could give large 
strong phases with QCD corrections, and hence large CP violation could be expected, so 
their effects cannot simply be neglected. However, in the QCDF method, the annihilation 
topologies (see Fig. 2) violate factorization because of the endpoint divergence. There is 
similar endpoint divergence when considering the chirally enhanced hard spectator scatter- 
ing. One possible way is to treat the endpoint divergence from different sources as different 



phenomenological parameters |22|. The corresponding price is the introduction of model 



dependence and extra numerical uncertainties in the decay amplitudes. In this work, we 



will follow the treatment of Ref. |^ and express the weak annihilation topological decay 
amplitudes as 

A-{B^PV)cxfBfpfvT.^ph , (2.5) 



where the parameters bi are collected in Sec. [II D|, and the expressions for the weak annihi 



lation decay amplitudes of B—>-PV are listed in the Appendix ^. To distinguish the decay 
amplitudes Eg. (p.4[) from Eq.( [2.5[ ), we will add a superscript / to the symbol in Eq. ([2.4[ ), 
and write it as A{B^PV)^A^ (B^PV). 



C. The QCD coefficients 

We now present the QCD coefficients in Eq.(p.4[). In our calculations, we neglect 
the mass of light quarks when we apply the equation of motion to the external quarks. 
We consider the chirally enhanced contributions from twist-3 LCDAs of the pseudoscalar 
mesons. As for the vector mesons, only the leading twist LCDAs for the longitudinally 
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polarized components $y contribute to the hadronic matrix elements, while the effects of 
twist-2 LCDAs for the transversely polarized components $y and higher twist LCDAs of 
the vector mesons are power suppressed and can therefore be neglected within the QCDF 
framework. We express the coefficients Oj in two parts, i.e., Oj = aij + aijj. The first term 
tti^i contains the "nonfactorizable" effects which are described by Fig.l(a-f), while the second 
part Oj // corresponds to the hard spectator scattering diagrams Fig.l(g-h). 

There are two different cases according to the final states. Case I is that the recoiled 
meson Mi is a vector meson, and the emitted meson M2 corresponds to a pseudoscalar 
meson, and vice versa for case II. For case I, we sum up the results for as follows: 



C21 



1 + 



as,/ 



1 + 



4tt 
Cpas 
in 

CpOis 



C3 + — [1 + 



Vm 
Vm 



nP - r 



^5,/ — C's 



9i 



Cpa^ 

J- H -. VM 

i - V r 

47r 



+ 



M 



6- 



An 



+ 



a?,! = C7 + 
^8,/ = C's + 

0-9,1 = Cg + 
^10,/ = Cio + 

where Cp 



Cs 
Nc 
CV 
Nc 
Cio 

Cg 
Nr 



1 - 



■ y r 



1-6 
1 - 

1 



An 

CpOts 



M 



An 

An 

CpOi, 



An 



+ 
Vm 
-Vm 



An Nr 



PiP,ew 
- M,3 

9n Nr '' 



+ 



■r,P,ew 

a Pm,2 
97r Nr ' 



01, // = 

02, // = 

03, // = 

04, // = 

05, // = 
oe,// = 
07,// = 

Os,// = 
Og,// = 

OlO,// 



nCpoii 

c 

nCpdi 

c 

nCpoii 

c 

nCpOis 
~Nf 

nC pois 

w 



■C2H{BMi,M2), 
■CiH{BMi,M2), 
■CiH{BMi,M2), 

CeH'iBMi,M2) 



nC pois 



CsH'{BMi,M2) 



nC poii 

c 

nCpa 



CioH{BMi,M2) 
CgH{BMi,M2] 



(2.6a) 
(2.6b) 
(2.6c) 
(2.6d) 
(2.6e) 

(2.6f) 
(2.6g) 
(2.6h) 
(2.6i) 

(2.6j) 



2Nc 



, and Nc = 3. The vertex parameters Vm and V^j result from Fig.l(a-d); 



the QCD penguin parameters ^ and the electroweak penguin parameters Pm'^1" result 
from Fig.l(e-f). The expressions for the penguin parameters are a little different from those 
in Here we consider the corrections of the electroweak penguins to Pmi, and the 



contributions of the QCD penguins to PmT- They can be written as 



pp 



pP 



M,2 



ryp,ew 
^M,3 



a 



|ln^ + 1 - Gm{0) - Gm(1) 

6 jJ, 6 



q=q' 

-2C,. 



fix 



(2.7a) 



Ci[l\n^ + l-GMis,)] + {Cs-la 
3 3 



|ln^ + ^-GM(0)-GAf(l) 
3 /i 3 



-In GrM(Sg) 



-2C, 



{C, + N^C,)[^\n'^ + ^-GM{s, 



r4 mfe 2 1^ 1^ /.N 



Q /I 

q=q' Z 6 fi 



1 — a; 



gin- + - - GMis,) 



q=q' 



4 mfc 
-111 — 
3 /i 



GmIsq 



K.G 



c0 7^. 



(2.7b) 



(2.7c) 



(2.7d) 



The definitions of quantities V^f, V^^, GM{sq), and GM{sq) can be found in and 
Sq = my ml. In the expressions for Pf^j and P^f^l", q' will run over all the active quarks at 
the scale /i = 0{mb), i.e., g' = u,d,s,c,b. The parameters H{BMi, M2) and H'{BMi, M2) 
in aj^//, which originate from hard gluon exchanges between the spectator quark and the 
emitted meson M2, are written as 

fsfv 



H{BV, P) 
H'{BV, P) 



™2 aB-^V 



fBf^ 



V 



™2 aB^V 



J J ^ X y 
J J t, X y 



(2.8a) 
(2.8b) 



For case II, except for the parameters of H{BMi, M2) and H'{BMi, M2), the expressions 
for tti are similar to those in case I. In particular, we would like to point out that because 
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{V\{qq)s±p\0) = 0, the contributions of the effective operators Oe.s to the hadronic matrix 
elements vanish, i.e., the terms that are related to a^ g disappear from the decay amplitudes 
for case II. As to the parameters H{BMi, M2) and H'{BMi, M2) in Oj,//, they are defined 
as 



H{BP, V) 
H'{BP, V) 



fsff 



fsfp 



dx I dy- 

J J t, 



X 



y 



+ 



m X y 



d^ / dx 
JO Jo 



^^y^B{^)^v{x)^^p{y) , 2iipx^l{y) 



X 



y 



m X y 



(2.9a) 
(2.9b) 



Now we would like to make some comments on 



1. It was shown in p3 that at leading order approximation are renormalization scale- 
independent, i.e., 

4^ = (^^6,8) and 1^ = (2 = 6,8) (2.10) 
d In/i d mfi 

where is the chirally enhanced factor, e.g., r^(/i) = mt,(^)[m (/^)] ^"^^ decay 

B^^K^ . It is formally 0{KQCD/''^h) power suppressed, but numerically large. The 
deviation of the data in Table H from Eq. ( |2.10|) can be reduced further when the higher 
order radiative corrections to the hadronic matrix elements are included. So in the 
following calculations, we will fix = m^. In addition, it has been proved that are 
free of gauge dependence. These two points guarantee that the decay amplitudes are 
physical. 

2. From Eq.( p^.6| ) we can see that, in the heavy quark limit, when the corrections at the 
order of as and a are neglected, we return to the NF approximation. 

3. "Nonfactorizable" effects appear at the order of or /and KqcD/^b- This reflects the 
fact that the final state interactions are either dominated by hard processes or power 
suppressed. Therefore the strong phases for the class-I (ai dominant) decays are 
small and furthermore perturbatively calculable in the heavy quark limit. Moreover, 
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from the data in Table ||, we find that the "nonfactorizable" effects contribute a 
large imaginary part to 02,7, i.e., the class-II (02 dominant) decays, especially the CP 
asymmetries, might be sensitive to the "nonfactorizable" corrections. 



D. The annihilation parameters bi 



The parameters of bi in Eq. (|2.5|) correspond to weak annihilation contributions. Now we 



give their expressions, which are analogous to those in p2 
Cf 



6i(Mi,M2) = -^CiAl(Mi,M2) 

c 

C. 



&2(Mi,M2) 



C2Al(Mi,M2), 



hiM,,M,) = -^{C3Al(Mi,M2) + C5A^(Mi,M2) + [C, + N,Ce\4iM^, M^)} 



iV2 



64 (Ml , M2) = { C^A\ (Ml , M2) + Ce A'2 (^1 , ^2 ) } , 



6r(Mi,M2) 



iV2 

Cf 

c 

C 



(2.11a) 
(2.11b) 
(2.11c) 



{C,A\{M^,M2) + CrAl{M,,M2)+[C^ + N,Cs]Ai{M,,M2)}, (2.11d) 



6r (Ml , M2) = ^ {Cio (Ml , M2) + Cs (Ml , M2) } . 



(2. lie) 



Here the current-current annihilation parameters 61^2 (Mi, M2) arise from the hadronic ma- 
trix elements of the effective operators Oi,2, the QCD penguin annihilation parameters 
&3,4(Mi,M2) from O3-Q, and the electroweak penguin annihilation parameters 63'^(Mi,M2) 
from Oj-iQ. The parameters of bi are closely related to the final states; they can also be 
divided into two different cases according to the final states. Case I is that Mi is a vector 
meson, and M2 is a pseudoscalar meson (here Mi and M2 are tagged in Fig. 2). Case II is 
that Ml corresponds to a pseudoscalar meson, and M2 corresponds to a vector meson. For 
case I, the definitions of y4^'-^(Mi, M2) in Eq. ( p.ll| ) are 



A(iV,P)=0, AiiV,P) = 0, 

AiiV, P) = na, f'dx f dy^y{x)¥'p{y)'^^'' + 
Jo Jo 



m x^y 



A\{V, P) = na, f'dx f' dy^v{x)<^p{y) 
Jo Jo 



1 1 

+ 



y{l- xy) x'^y 

1 /■! r 1 1 

A\{V,P) = -Tias I dx / dy^v{x)^p{y) 



x{l — xy) xy"^ 



(2.12a) 
(2.12b) 

(2.12c) 

(2.12d) 
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Al{V,P) 
For case- II, 

A({P,V) 

4{p,v) 

A\{P,V) 
A^{P,V) 
AliP,V) 

Here our notation and convention are the same as those in |2^. The superscripts i and / 
on A''-^ correspond to the contributions from Fig.2(a-b) and Fig.2(c-d), respectively. The 
subscripts A; = 1, 2, 3 on A*^'-^ refer to the Dirac structures (V — A)®{V — A), {V — A)(Si{V + A) 
and {—2){S — P)^{S + P), respectively. ^v{x) denotes the leading-twist LCDAs of a vector 
meson, and ^p{x) and $p(a;) denote twist-2 and twist-3 LCDAs of a pseudoscalar meson, 
respectively^ 

Note that assuming SU{3) flavor symmetry and symmetric (under x^x) LCDAs of light 
mesons, we have A\ = —A\. In this approximation the weak annihilation contributions (for 
case I) can be parametrized as 



Tcasf dxf dy^v{xWp{y) 
Jo Jo 



m xy{l — xy) 



(2.12e) 



0, 4{P,V) = 0, 

-na, I'dx /'«(a:)$v(l/)^^^^: 
J J m XV 



nag dx dy^p{x)^viy) 
Jo Jo 



+ 



y{l - xy) x'^y 

1 /•! r 1 1 ] 

Tcas I dx dy^p{x)^v{y) -r. ^ + — T 

JO ^x[l — xy) xy^^ 

m xy{l — xy) 



nas dx dy(^^p{x)(^viy) 
Jo Jo 



(2.13a) 
(2.13b) 

(2.13c) 

(2.13d) 

(2.13e) 



^ The expression for is different from that in |24|. This difference originates from the two 
different ways of dealing with the annihilation contribution from the twist-3 LCDAs. Cheng and 
Yang calculated them in coordinate space, and we take the projection in momentum space given 
in |2^. We made some investigation of this problem in |27], and discussed it with Beneke et al. 
The point lies in how to deal with the surface terms in the integral. In addition, we found that, 
when the annihilation contributions from the twist-3 LCDAs are considered for B^PP decays in 
coordinate space, not only logarithmic divergence but also linear and quadratic divergences appear, 
which is a very serious problem. So we use the method given in |2^, which gives only a logarithmic 
divergence even for B^PP decays. 
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(2.14b) 



27r^-6(Xl + 2X^^ 
4(\/,P)~67ra,r^(2Xl-XA), (2.14c) 

where = J^dx/x parametrizes the divergent endpoint integrals. We can get similar 
forms to Eg. ( p. 141) for case II, but with Al{P,V) = -Al{V,P). In our calculation, we will 
treat as a phenomenological parameter, and take the same value for all annihilation 
terms, although this approximation is crude and there is no known physical argument for 
justifying this assumption. We shall see below that Xa gives large uncertainties in the 
theoretical prediction. 



III. INPUT PARAMETERS 

The QCDF expressions for the hadronic matrix elements are written as the product of 
hard-scattering kernels and LCDAs of mesons. The hard scattering kernels are perturba- 
tively calculable, while the soft and nonperturbative effects are incorporated into two kinds 
of universal parameter: the form factors and the LCDAs of mesons. The decay amplitudes 
are also related to the CKM matrix elements, the masses of quarks and mesons, various 
decay constants of mesons, and so on. We will specify them in the following discussions. 



A. The CKM matrix elements 



The CKM matrix is described by four independent parameters. Using the Wolfenstein 



parametrization ||3^, it can be expressed as 

/ 



VcKM 



AV2 A A\\p - if]) 

-A 1 - AV2 AA2 

_ ^ _ — AA^ 1 



^ y4A^(l - p-iT]) 



C(A^). 



(3.1) 



Two of them are well-determined: A = 0.819±0.040 and A = 0.2237±0.0033 PB|. Recently 



M. Ciuchini revised the values of p and rj in p = 0.218±0.038, rj = 0.316±0.040, and 
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7 = (55.5±6.2)°. This gives p = 0.224±0.039 and r] = 0.324±0.039. A. Hocker, et al. |38 
advocate another approach to a global CKM matrix analysis, using the frequentist statistics 
named RRt; their best fit results are A = 0.83±0.07, A = 0.222±0.004, p = 0.21±0.12, 
fj = 0.38±0.11 and 7 = (62±15)° at 95% confidence level. H. Lacker and M. Neubert et al. 
suggest that it is possible to derive constraints on 7 from a global i?fit analysis of various 
CP-averaged branching ratios B^mr, ttK; their best fit of the QCDF theory to the data 
yields {p,f]) = (0.05,0.381) with x^/udof = 0.46 [^, and correspondingly 7^83°. In our 



calculation, we will take the central values of the results in as input parameters. At the 



same time, we will also give estimations based on the results in 39 



B. Quark masses 

The masses of quarks appear in our calculations in two different ways. The pole masses of 
quarks arise from the loop integration over the virtual internal quarks in penguin correction 
diagrams. They contribute to the penguin parameters and Plf^ in Eq.([^.7D in terms 
of GM{sq) and GM{sq) where Sq = niq/ml. We fix them as 



f^u = = = 0, = 1.45GeV, nih = 4.6GeV. (3.2) 

The other type of quark mass is the running quark mass which is renormalization scale- 
dependent. It appears in terms of the chirally enhanced factor which arises from the 
hadronic matrix elements of {S + P)^{S — P) operators through the equations of motion, 
and the twist-3 LCDAs of the pseudoscalar mesons. Estimations of the running quark masses 
within the QCD sum rules approach are collected in |40]. Here we would like to use the 



central values of Particle Data Group 2000 data gT| for discussion. 

mbiffib) = 4.2GeV, rn,(2GeV) = 122.5MeV, (3.3a) 

md(2GeV) = 6MeV, m„(2GeV) = 3MeV. (3.3b) 

Using the renormalization group equation, we can get their corresponding values at the 
scale p = 0{mfy). In addition, we would like to point out that, because the running masses 
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of light quarks have large uncertainties, we will take the r'^'\l — /^(,)//^(,))~r^^r^=r^ 
approximation for simplicity in our numerical calculations. 

C. The form factors and decay constants 

Now let us parametrize the hadronic matrix elements of {PV\Oi\B) j in Eq. (p.4|) . After 
Fierz reordering, with the NF assumption |jl8|, they can be written as 



{PV\0,\B) = Z,{P\r\0){V\J,\B) + Z2{V\J'^\0){P\J'^\B), (3.4) 

where and J''^ are hadronic currents and Zi 2 are the corresponding coefficients. The 
hadronic current matrix elements are defined as follows [p^ : 

{Pmri^q\0) = -ifpl^ {V{l,e)\q^^q\Q) = /ymye*^ (3.5a) 

{Pmra - 75)g|S) = K + ^'- ^^~^^ fc^]Fi(A:^) + "^^~"^H ^Fo(A;^), (3.5b) 

rriB + my k'^ 
e*-k 



+ ie*Jmb + mv)Ai{k^) - i {pb + l)^,A2ik'') 

^ mb + my 

2niv{e*-k) 



„ ^2 -k,Mn. (3.5c) 

where k = pb — I, and e* denotes the polarization vector of the vector meson V. fp and fv 
are decay constants, and Fo,i(A;^), K(fc^), and Ao,i,2,3(fc^) are form factors. In addition, at 
the poles /c^ = 0, we have 

Fo(0) = Fi(0), A(0)=A3(0), (3.6a) 

2myA3(0) = {niB + my)Ai(O) - (m^ - 7/1^)^2(0). (3.6b) 

In the QCDF framework, they are all nonperturbative quantities, and appear as universal 
input parameters. As a good approximation, we take these form factors at /c^ = in our 
calculations. 

Because the fiavor octet mixes with the fiavor singlet in the SU (3) quark representation 
of mesons, the corresponding hadronic parameters are difficult to determine. Here we assume 
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ideal mixing between oj and 0, i.e., oj = {uu + dd) / \/2 and (j) = ss. As to r] and r]', we follow 
the convention in p9|j4^P3[] : the two- mixing- angle formula is applied for the decay constants 
and form factors, and the charm quark content in rj and rj' is assumed to be negligible: 

{0\qi^l5q\v'-'\p)) = {q = u, d, s) (3.7a) 



2 



<" = f£"(^-5ij). F-;' = F„-!C^^^). (3.7e) 

We take the mixing angles as 6^ = —22.2° and = —9.1°. For the B^K transition form 
factors we use the SU{'i) flavor symmetry approximation 

The decay constants and form factors are nonperturbative parameters; they are available 
from the experimental data and/or estimated with well-founded theories, such as lattice cal- 



culations, QCD sum rules etc. Now we take their values in our calculations as [B3,0,|l-H^: 



U = 131MeV, fK = leOMeV, fx- = 214MeV, 
fp = 210MeV, = 195MeV, = 233MeV, 

/o = 157MeV, fs = 168MeV, fe = ISOMeV, 
F^J{0) = 0.28±0.05, ^^^*(0) = 0.39±0.10, 

Ao''{0) = 0.30±0.05, ^^""(0) = 0.30±0.05. 



D. The LCDAs of mesons 



The LCDAs of mesons are also basic input parameters in the QCDF formula Eq.( |2.3| ). 
The LCDAs of light pseudoscalar mesons are defined as |]48[| : 
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{PikMz2)qiz,)\0) 



(3.9) 



(here mi and 



where fp is a decay constant, z = Z2 — Zi, x = 1 — x, and /ip = yrnx{ti)+m2{iJL)] 
7712 are running masses of the valence quarks of the pseudoscalar meson (P) with mass mp.), 
$p(a;) is the leading twist LCDAs, and $p(a;), $p(x) are twist-3 LCDAs of the meson. We 
shall use the asymptotic forms of the LCDAs for the following discussion: 



<l>p(x) = Qxx, $^(x) = 1, ^pix) = Qxx. 



(3.10) 



The twist-2 LCDAs of the vector mesons are defined as 47,49 : 



{0\m^,Mz)\Vik, A)) = tie^^K - etk,)f^ / ^dx e 



-ixk-z ^1. 



{0\miM^)\Vik,\)) 



k^^fvmv f dxe "■^"<l>{.(x) 
k-z J 



(3.11a) 
(3.11b) 



where e is a polarization vector, and for longitudinally polarized mesons ey = k/ray. ^y{x) 
and ^y{x) describe the quark distributions of transversely and longitudinally polarized 
mesons. In our calculations, the contributions from ^y{x) are power suppressed and hence 
can be neglected, i.e., we take the following approximation: 



^y{x) = ^y{x) = Qxx. 



(3.12) 



For the wave function of the B meson, we take the form in p^ , |35| , |50[1 : 



$B(0 = A^Br(l-0 exp 



2ujI 



(3.13) 



where A^p is the normalization constant. $p(0 is peaked around ^~0.1 with uob = 0.4GeV. 

With these LCDAs of the mesons, we still cannot get the parameters and hi at once, be- 
cause we will encounter endpoint divergence in dealing with the integrals Eq.( p.9| ), Eq. ( p.l2| ), 
and Eq.( p.l3D . If the transverse momentum k_[_ of the partons and the Sudakov resumma- 
tion can be taken into account consistently in the QCDF approach, the above-mentioned 
integrals might be convergent. But now we have to treat these divergent integrals as phe- 
nomenological parameters: 



where q can vary from to 6, is an arbitrary phase, and O°<0<36O°. These notations are 
almost the same as those in |2^j2^ . In our calculation, and Xa denote the endpoint di- 



vergent integrals from hard spectator scattering Eq. (|2.9| ), and weak annihilations Eg. ( p.l2| ), 
Eq.( p.l3|) , respectively. We take A = Aqcd and ge~^'^ = in as default values. 

IV. BRANCHING RATIOS 

The branching ratios of charmless decays B-^PV in the B meson rest frame can be 
written as 

Br(B^PV) = —^\A(B^PV)\\ (4.1) 
In our calculation, we take r^o = 1.548ps, rp± = 1.653ps, and 



J[m% - {mp + my)2][m| - (mp - mvY^ ^ 
1^1 = 2^, • ^'-2) 

Since QCD factorization works in the heavy quark limit, the above masses of light mesons 
should be taken as zero for consistency. 



The experimental measurements are collected in Table |T|. Our numerical results for CP- 
averaged branching ratios for B^PV are listed in Table |^ and Table 0, calculated at the 
scale /i = mf, for two choices of CKM matrix elements: (1) A = 0.819, A = 0.2237, p = 0.218, 



and f] = 0.316 0; (2) A = 0.83, A = 0.222, p = 0.05, and r/ = 0.381 H. The other related 



parameters are taken at their default values. In the tables, Br stands for the values that 
are calculated without including weak annihilation contributions within the NF framework; 
Br-^ (x\A^ and Br^^"'(x\A^ ^A^]"^ are estimated with the QCDF approach where A^ denotes 
the annihilation contributions. Since the above two choices of CKM matrix elements are 
nearly equivalent to two choices of the angle 7, many of the CP-averaged branching ratios 
are similar to each other for these two sets of CKM matrix elements, as can be seen from 
Table |V| and Table 0. Only the decay channels that have large interference between tree 
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and penguin contributions are sensitive to the angle 7, such as ^tt^ K*~ and -^-k^ K*~ . 
It is thus possible to extract 7 from these decay channels. However, we must consider the 
uncertainties due to the variations of various parameters, such as form factors, and the 
model dependence of the parametrization of chirally enhanced hard-spectator contributions 
and annihilation contributions, and so on. Nevertheless, it is still possible to obtain some 
information on 7 with detailed analysis. The readers may notice that some numerical results 
in Table ^ and Table |V] are inconsistent with the experimental measurements listed in Table 
[m| . The inconsistency is not serious because we use default values of the parameters for 
the tables and do not consider the uncertainties from the input parameters. In fact, we will 
see that for appropriate regions of parameters predictions from the QCDF approach are in 
agreement with current measurements for most of the B^PV decays. 

From the experience m B PP analysis |2^, we know that the main theoretical 
uncertainties of branching ratios come from CKM matrix elements, form factors, and weak 
annihilation contributions. It is easy to imagine the importance of CKM matrix elements and 
form factors, but for weak annihilation topologies the importance was first noticed recently 



m 



34 1 within the perturbative QCD approach. However, within the QCDF framework, the 



annihilation topologies introduce endpoint divergence, which violates factorization; in Ref. 



22 1 , the authors phenomenologically parametrize the endpoint divergence integral as Xa- 



In this work, we will follow their way of estimating the annihilation contributions, although 
it will introduce model dependence and numerical uncertainties. 

We know that the leading power of the annihilation contribution is always X\ for B^PP 
decay amplitudes. But for B^PV, the case is different: According to the formulas in Sec. 
|II D| , b^^ and bf^ are negligible in general because of small Wilson coefficients, and 64 is also 
small due to the cancellation between A\ and Ag. For the other bi parameters we have 

b, = ^C,A\ = 18^C,na,{XA + y - 4), (4.3) 
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h ~ ^N.CeM = ±Q^N,Cer^7casi2Xl - Xa). (4.5) 

So when bi or 62 dominates the leading power of the annihilation contributions is Xa- 
This indicates that in some cases, compared with B PP decays (where appears), 
annihilation topologies might introduce smaller uncertainties into B^PV branching ratios. 
In the following we will proceed to analyze some decay modes for B^PV in detail. 

A. B'^^'k^p^ decays 



From Table |ITI| and Table |V|, we can see that the theoretical estimation of the CP- 
averaged branching ratio Br{B^^7i^p^ + 7r^p+) is in agreement with the measurements of 
the BaBar and CLEO Collaborations within one standard deviation, especially when we 
consider the uncertainties due to the variations of input parameters (see Fig.^. 

The decays of B^^-K^p^ are ai dominant. The "nonfactorizable" contributions are small, 
so there is no distinct difference between the results obtained with the QCDF approach and 
those with the NF approach. We can see this point from Table fV[ The annihilation ampli- 
tudes of these decay channels are hi dominant, so they do not contribute large uncertainties 
to the branching ratios. From Fig. ^(c), it is obvious that the hard spectator scattering has 
little impact on the branching ratio. Therefore the main uncertainties in the CP-averaged 
branching ratios of these decay modes originate from the form factors and CKM matrix 
elements (especially |V^„fe|), as can be seen from Fig.|^. 

The ratio of branching ratios g^^^oZ^^+p-j has been discussed in [|9|. From Fig.|](f), 



it is obvious that Br{B^^7i^ p^) > Br{B^^7c^p~), mainly due to the difference between 
the decay constants: fp > f-,^. In addition, the destructive interference among penguin 
amplitudes also lowers the branching ratio of B^^n^p^, although this effect is relatively 
small. In fact, because the decay constants of vector mesons are generally larger than those 
of pseudoscalar mesons, i.e., fv > fp, it seems universal that, for the ai dominant B^PV 
decay modes, the decay channels with emitted vector mesons have larger branching ratios 
than the corresponding decay channels with emitted pseudoscalar mesons. From Fig.^(a), 
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it is clear that the averaged branching ratio of these decays is only mildly dependent on the 
angle 7, which hints the small contribution from penguin amplitudes. So it is reasonable to 
neglect the penguin amplitudes, which leads to 



(4.6) 



Clearly, this ratio is insensitive to the CKM matrix elements, dynamical coefficients Oj, and 
so on. Since it relates different form factors, a measurement of this ratio may be helpful to 
improve theoretical predictability. We draw this ratio versus 7 in Fig.|(f). 



B. B — >7rp, TTu; decays 

The decays B~ ^Tip^nuo are also tree dominant, but they are determined by Oi + ^02. 
From Fig.^(a) and Fig.|^(a), we can see that within appropriate ranges of the parameters, our 
results for Br{B^ ^ti^ p^) and Br{B^ ^ti^uj) are in good agreement with the measurements. 
For BrlB^^n^p"), the branching ratio is predicted numerically to be large and it should 
be observed readily, although there is only an experimental upper limit so far. 

Generally speaking, for the branching ratios, the numerical uncertainties due to the 
variations of the CKM matrix elements and form factors are always large. But it is distinctive 
for B^^np, -Kuj decays that hard spectator scattering also causes sizable uncertainties [about 
20%; see Fig. H(c) and Fig. ||(c)]. This is because, as mentioned above, these decay 
modes are determined by ai + and "nonfactorizable" effects; in particular, terms of 
hard spectator scattering contribute greatly to 02. As for annihilation contributions, due 
to the cancellations such as h2{P, V) - 62(V, P) and h{P, V) + hiy, P), A" {p- ^ti- p^) and 
A"- {B~ ^Ti~ uo) are dominated by 63 and 62, respectively. So for B~^7i~p^ the annihilation 
topologies also contribute large uncertainties to the branching ratio, while for B^^ti^uj, the 
uncertainty due to the annihilation parameter Xa is negligible [see Fig ^(b) and Fig|^(b)]. 



22 



C. B^K(f) decays 

The decays B^^K^ip and B~-^K~(f) have triggered great theoretical interest 
2^ , |29|j30| , |5T|| because they are pure penguin processes. Experimentally, these two decay 



modes have been observed by the BaBar, Belle and CLEO Collaborations; the averaged mea- 
surements ignoring correlations are Br{B^~->-K^(j)) = (7.5±1.8)x 10^^ and Br{B^^K^(j)) = 
(7.9±1.2) X 10^^. Our results within the QCDF framework are in good agreement with the 
experimental data within la [see Fig.^(a,b)]. 

From the expressions for the decay amplitudes in Appendix B of , we can see that the 
decays B^Kcf) only relate to the CKM factor Vi^Vj*, which is well determined experimentally. 
So the decay rates are independent of the angle 7 and the uncertainties due to the variations 
of the CKM parameters should be very small. 

For B^K(f) decays, it is very interesting that the prediction is much different between 
the generalized factorization (GF) approach and the QCDF approach. Within the GF 
framework, the predicted branching ratios of the decays B-^Kcf) are very sensitive to the 
"nonfactorizable" effects, they vary from 18x10"^ to 0.4x10"^ p9[ (or from 13x10"® to 
0.3x10"*' 0) with the effective number of colors N^^f = 2~oo. But within the QCDF 



framework, there is no need to introduce the effective color number N^^^ . However, since 
A"'{B^K(f)) is dominated by 63, the branching ratios will be very sensitive to the annihilation 
parameter Xa- As for hard spectator scattering, since {B^Kcj)) is dominated by 04, the 
uncertainties from the variations of Xh should be small. In ||2^, the authors also analyzed 



the decays B^Kip with the QCDF approach. The effects of the twist-3 LCDAs of K mesons 
on the hard spectator interactions were considered there in spite of their small corrections 
to the decay amplitudes. When the contributions from the weak annihilations are taken 
into account, their predictions are Br{B'^^K^(j)) = (4.0ti;4) x 10"® ^-^d Br{B^^K^(j)) = 
(4.3li;4) X 10"® (their central values correspond to Qh = Qa = 0), which are consistent with 
our estimation at = m?,, although the values of input parameters are slightly different. 
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D. B^rjK* decays 

The B^rjK* decays are well established experimentally by the BaBar, Belle, and CLEO 
Collaborations (see Table [m| ). But within the QCDF framework our estimations (see Table 
0) using default input parameters are several times smaller than the experimental measure- 
ments. Even considering the large numerical uncertainties, the predicted branching ratios 
are only marginally consistent with the experiments [see Fig. |^(c,d)]. 

In fact, there is a similar problem for B^rj^'^K decays which aroused intense discussion 
several years ago (see, for example, Refs. [^^-|57|). It is now commonly believed that this 
puzzle may be related to a special property of t]': it has large coupling with two gluons 
Il55| , |56|| . However, there are definitely no large rigg contributions to the decay amplitudes of 

Here we would like to point out that the CP-averaged branching ratios for B—>-r]'K* 
decays are similar to those for B-^rjK* decays with the default parameters (see Table 0), 
which are different from the results listed in p0| , p9| , |57| . H. Cheng and K. Yang thought that 



Br{B—>-r]K*) should be much larger than Br{B—>-ri'K*) in the GF approach, and their favor- 
able values are ^^g^|S^~20 and ^|gl^^a^^l8 with N^f^iLL) = 2 and N'J^{LR) = 6 



|20|. In [E^, A. Ah et al. predicted that Br{B^r]'K*) would be too small to exceed 1x10" 



In those works [^,^,0], the small branching ratios for B^rj'K* decays are due to delicate 
cancellation between different parts of the decay amplitudes. However, this cancellationis 
sensitive to the choice of the input parameters such asform factors and so on (see Table. IV 
in [^). In our analyses, we find that the weak annihilation contributions might cause large 
uncertainties in the predictions of B^r]^'^K* decays. 

Let us have a closer look at the numerical uncertainties of Br{B—>-riK*). Their decay 
amplitudes are penguin dominant, so for the CKM factors only iVf^V^*! plays an important 
role. Since |Vtf,VJ*| is well determined experimentally, we do not expect large uncertainties 
from the CKM parameters. As for hard spectator scattering, since 04 and are the most 
important coefficients in the amplitudes, the hard spectator parameter Xh has a small effect 
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on the branching ratios. The decay amphtudes also depend on the form factors Ff^^ and 
A^^^\ but it seems unhkely that the form factors would be dramatically different from the 
default values. Fortunately, the annihilation contribution is dominated by 63, which could 
bring large uncertainties to the branching ratios. However, considering the big gap between 
experiments and default numerical predictions, the annihilation parameter should be 
very large. In this paper, we choose the model that Xj^ is universal for all B—>-PV channels, 



and even universal for all B^PP channels |28], so is it possible for a large Xa to survive 



for a global B-^PV and PP fit? It may be a challenge for the QCDF approach to obtain a 
large Br{B^rjK*). 

E. Other decay modes 

There are still some decay modes for B^PV where the QCDF predictions are marginally 

— *o 

consistent or even inconsistent with the experimental results, such as B~^tx^K and 
B^^K' ,7i^K*^ (see Fig. H). As an illustration, let us see what the QCDF predic- 
tion for B^^n^K and its uncertainties are. Without considering annihilation topologies, 
this decay channel is a pure penguin process and dominated by the QCD coefficient and 
the CKM factor |V(feV^*|. So the branching ratio is independent of the angle 7, and the 
variation of the CKM parameters and hard spectator parameter will have very little effect 
on the decay rate. As for the form factors, only F^~*'^ is involved. But from the experience 
of B^TiTi ||28| we know that it is favorable to have a small form factor. Hence we need 



a large annihilation parameter X^ again. This indicates that soft interactions play a very 
important role in some cases. Anyway, it is not good news because the soft annihilation part 
is model dependent. As we know, annihilation topologies are very important for B^PP 
decays; therefore it is very important to check whether a large Xa could be acceptable for 
other decay modes, such as B^-k-k and nK. Experimentally, it is also important to search 
for pure annihilation processes, such as B^—>-K^K~ and K^K*^, which may be helpful in 
learning more about the annihilation mechanism. 
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Within the QCDF framework, some B^PV decay modes are predicted to have 
small branching ratios. One type is the 02 dominant neutral B decays, such as 
i?^^7r'^p°, TT^cij, r/'^'^p", r/'-'^u; etc. Another type is those decays whose tree amplitudes are 
denoted by a2 and suppressed by the CKM factors; moreover, the interference among pen- 
guin amplitudes is destructive, for example, the decays ^—yK^p^,TC^uj can be classified 
into this type. The third type is pure penguin processes including B^^l(^ , K^K*^ and 
K*^ ] their small branching ratios are due to the delicate cancellations among various 
competing terms. The last type is those decays that are governed by the small coefficients Oj, 
or electroweak penguin dominant decays, including B~^(j)Ti~ and 'B^ —^(pTT^ , (j)ri^'\ etc. Our 
numerical results in Table and Table indicate that the CP-averaged branching ratios 
of the above-mentioned decay modes are indeed small, and do not exceed lxlO~^ without 
considering the uncertainties. They are much below the corresponding experimental upper 
limits. But for these decays, the effects of weak annihilations, soft final-state interactions, 
and other kinds of power corrections might be very important or even dominant. 



V. CP VIOLATION 

CP violation and quark mixing are closely related to each other in the SM. The existence 
of a CP- violating phase in the top sector of the CKM matrix has been established, i.e., 
lm{Vtd)o:fj^O. Large CP-violating effects are anticipated and have been observed in the B 
meson system. In this section, we will study CP violation for charmless decays of B—^PV 
within the QCDF approach. 

Suppose that the decay amplitude for B^f can be expressed as 

A{B^f) = Aie-'^'e-'^' + A^e'"'^^ e-'^\ (5.1) 

where the Ai are the magnitudes, and and 5i are the strong and weak phases, respectively. 

Then we can get the CP-violating asymmetries 

_ V{B^f) - T(B-^f) _ lAjB^f)]' - \A(B-.f)\' 
T{B^f) + T{B^f) 1^(5^/) P + \A{B-.f)\^ 
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2AiA2sin(^i - 6)sin(5i - ^2) 2) 



AI + AI + 2AiA2Cos(6 - 6)cos((5i - ^2) ' 
clearly, ^c'pocsin(^i — .^2)- So we would like to repeat some remarks. (1) Within the NF 
framework, there are no direct CP violations for the hadronic charmless two-body B decays, 
because in the decay amplitudes the Wilson coefficients and the hadronic matrix elements 
which are expressed by the product of the decay constants and form factors, are all real, so 
the strong phase shift ^1— ^ = 0. This indicates that "nonfactorizable" effects are important 
for CP violations in hadronic B decays. (2) In the QCDF framework, the CP violations for 
the class I (ai dominant) decay modes should be small, because the strong phases arise at 
the order of and/or h^Qco/^b) and hence sin(^i — ^2) is small in general. For class II 
(02 dominant) decay modes, since 02,/ has a large imaginary part (see Table |I|), a large CP 
violation might occur. 

For charged B decays, the final states are self-tagging. The direct CP-violating asymme- 
tries are defined as 

- r(i?-^/-)-r(5+^/+)- ^^-^^ 

For neutral B mesons, effects of oscillation will induce CP asymmetries, so time- 

dependent measurements of CP-violating asymmetries are needed: 

^ ngW^/VrW)-/). (5.4) 
r(B°(i)^/) + r(i?»(«)^/) 



As discussed in |p9|,p8[, here we consider three cases for decays of B and B^^PV . 



case 1: 5°^/, W^f, but B^^f, Wy^f, for example, W^K'p+. CP-violating 
asymmetries for these decays are similar to those for i?^ decays, and no mixing is 
involved. 

case 2: = /)<— B^, for example, 'B^&z:B^^TT^p^,riuj, etc. The time-integrated 

asymmetries are 

1 X 

•^CP = T"^ — 2^^' + T"^ — ^0,+,', (5.5) 
1 + 1 + 
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1 — |Acpp 

i + iAcpr 



-2 Im(Acp) _l^^42_(0H/) 



1 + |AcpP 



(5.6) 



where x = Am/T = 0.73±0.03 a^/ and a^+e' are direct and mixing- induced 

CP-violating asymmetries, respectively. 

case 3: B^^{f and There are three examples for B—>-PV decays; they are 5^ 

and B°^n^p^,K^K*^ and B° and B^^K^K*^ ,K%K*^ . We follow the conventions 



for the time- dependent asymmetries in p9| , p8| . The four basic decay amplitudes of 
transitions and B^{t)^f and / at t = are written as 



g = A{B'{0)^f), 
h = A(B'{0)-.f), 



h = A{B\0)^f), 



(5.7a) 
(5.7b) 



For example, if / = K^K* , then h and g will correspond to Eg . (|A2|) and Eq. 
The time-dependent decay widths are written as 



T{B\t)- 


-/) = 


2 V'^ 


^+i/^r) 


1 + Og/ cos (Ami) + 


ae+e'sin(Amt) 




-/) = 


2 V'^ 




1 — %'Cos(Amt) — 


a^+^/sin(Amt) 


r(fi°(t)- 


-/) = 


2 




1 + a^'Cos(Amt) + 


ae4.e'sin(Amt) 


T(B^{t)- 


-/) = 


e-r* . 
2 i'^ 




1 — ae/cos(Amt) — 


ae+e/sin(Amt) 



where q/p 



and 



Of' 





\h\^ 








\9\' 


\h\^ + 


\9\' 



-2 lm[iq/p)x{h/g)] 

l + \h/g\' 
-2 lm[iq/p)x{g/h)] 

i + \g/h\' 



(5.8a) 
(5.8b) 
(5.8c) 
(5.8d) 



(5.9a) 
(5.9b) 



Our results for CP- violating asymmetries for decays B-^PV are listed in Tables |VJ-|VTT|. 
The parameters (including a^/, a^', a^+e', o-t+i'^ and Acp) with superscripts / and f + a are 
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calculated in the QCDF framework with decay amplitudes A-^{B^PV) and A-^ (B^PV) + 
A°'{B^PV), respectively. Some remarks are in order. 

(1) Within the QCDF framework, the strong phases are either at the order of a.^ or power 
suppressed in KqcD/'^nb- Only radiative corrections are in principle calculable in the QCDF 
approach. However, numerically power corrections in KqcD/^b might be as important as 
the radiative corrections. So in a sense the QCDF method can predict only the order of 
magnitude of the CP-violating asymmetries. 

(2) The theoretical predictions for CP-violating asymmetries are compatible with mea- 
surements within one standard deviation (see Fig. However, so far, the present 
experimental measurements on B^PV have too large uncertainties, e.g., Acp{'^t^^) = 
-0.34±0.25 or -0.0li[J:^?±0.03 0, and Acp{<pK^) = -0.05±0.20±0.03 ig, so they 
cannot provide any useful information. 

(3) Although the uncertainties from variations of nonperturbative parameters, such as 
form factors and so on, are reduced to some extent for the CP-violating asymmetries, the 
weak annihilations can still have great effects on the CP asymmetries for some decay chan- 
nels, such as a^i and a^^^i for the 02 dominant decays 'B^^nu, rjp^, etc. 

(4) It is worth noting that for the decays B'^^cpir, (f)r]^'\ there are no direct CP asym- 
metries a^i within the QCDF approach, while in the GF framework, a^' varies from ~1% to 
-23% with N^ff = 2~oo |29|. 

(5) One might wonder why the CP- violating asymmetries for the Oi dominant decays 
^o_j^^±p=F gj-Q not as small as we expected. In fact, we only expect that the direct CP 
asymmetries at t = are small, i.e., a^// = and a^" = jj^p'jllp . For instance, we 
take g = A{B^^7i- p^); then with the CKM parameters in af„ = 2.3%, af,;^" = 5.2%, 
at, = —0.2% and a£^" = 4.6% at the scale /i = m^. They are really small. 
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VI. CONCLUSION 



(1) In the heavy quark hmit, we calculated the hadronic charmless decays B^PV with 
the QCDF approach including chirally enhanced corrections and weak annihilation topolo- 
gies. Neglecting the power suppressed AgcD/'mb effects, the "nonfactorizable" contributions 
that cannot be calculated within the NF framework are perturbatively calculable with the 
QCDF approach at least at the order of ccg. They provide the strong phases which are 
important for CP violations. 

(2) Most of the CP-averaged branching ratios are in agreement with present measure- 
ments, but there are very large theoretical uncertainties due to the variations of input 
parameters. However, some decay channels, such as B^rjK* and B~^7r~K , are only 
marginally consistent with the experimental observations. In these decay channels, nonper- 
turbative contributions, such as weak annihilation topologies, play a crucial role and need 
further investigation. 

(3) The CP asymmetries for the oi dominant decays are small because of small strong 
phases, 0{as) and/or 0{!s.QCD/'mb)- However, the "nonfactorizable" effects and the weak 
annihilations can provide large imaginary parts for the 02 dominant decay amplitudes, and 
lead to large direct CP asymmetries. The theoretical predictions for CP-violating asymme- 
tries are compatible with current measurements within one standard deviation. It is worth 
noting that the QCDF predictions might give only the proper order of magnitude of the CP 
asymmetries. 

(4) The present experimental data for nonleptonic B meson decays are not yet sufficient, 
especially for measurements of CP violation. On the other hand, there are still many un- 
certainties in the theoretical frame, for instance, the weak annihilations and other potential 
power corrections. Great advances in both experiment and theory in the near future are 
strongly expected. 
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APPENDIX A: THE ANNIHILATION AMPLITUDES FOR B^PV 



1 

2 



2 



-*0 J^o^ 



(Al) 



V2 

HVubV:, + KbKd) [h{K+, K*-) + b,{K*-, K+) 
+br{K+, K*-) - ^br{K*-, K+)] |, 



(A2) 



V2 



(A3) 



V2 

-\bT{p-,T^^) + hT{T^^,p-) - \hT{p-y 



(A4) 



-lbr{K*-,n^)]y 



(A5) 
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*+ 



+bT{K*+,K-) - ^br{K-,K*+) 



V2 



fBfKfK*[iV, 



2 2 ~ 



a^b'^k'p') = -^fBfKf,{{v.bv:, + v,bV:,)[b,{K' , p') - hr(if,p'] 



A'^iB'^n'p') = §^fBUU{ Vu,V:U blip', + p') 



+2b,{n', p') + 2b,{p\ n') - ^b^p', 7r°) 



A'iB'^n'u;) = ^fj,fJ^!^V^,v:,[b,iu;,7r') + b,i7r',u;) 

HVubV:, + KbKd) [ - &3(7r°, a;) - bs{uj, 7r°) 
+lbr{n',u;) + hr{u;,n') + lbr{n',u;) 

+lbr{^y) 
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^"(S^^ttV) = 0, 



(17*^ 



+2h(r)^'\u;) + 2h(u;,rj^'^)-lbr(r)^'\u;) 



*0 



+(K6K:<i + Wc^) lh(K', K*-) + br{K', K*-)]y 
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A^B-^TT-cu) = ^fBfnU[VubV:,[b2{n-,u;) + b2{u;,n-) 
+6r(7r-,cu) + 6r(^,vr-)" 



A\B-^t:-4)) = 0, 



+6r(p-,vr°)-6r(7r°,p-) 



+hT{p-.V^'^)-rbT{v^'\p- 



A\B-^r]^'^K*-) = 



+ KbK^) bs{K*", K-) + 6r (i^*", i^-) 
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V2 

+{Vubv:, + k^k:) \h<i>: K-) + 6^(0, K-) 



(A30) 



A\B-^K-u) = ^fBfKUV^bV:MK-,u;) 



(A31) 



(A32) 



^/B/i./p{v;602(^",P") 



HVu,V:, + K^K:) p-) + hT{K\ p-)] }. (A33) 
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FIGURES 




(b) (c) (d) 




(f) (g) (h) 



FIG. 1. Order of corrections to the hard-scattering kernels. The upward quark hnes repre- 
sent the emitted mesons from b quark weak decays. These diagrams are commonly called vertex 
corrections, penguin corrections, and hard spectator scattering diagrams for Fig.(a-d), Fig.(e-f), 
and Fig.(g-h) respectively. 
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FIG. 3. Decays W and -^TT^p^ versus 7 at the scale = mi, within the QCDF approach, 
including the effects of weak annihilations. Fig.(a-e) are CP-averaged branching ratios, and Fig.(f) 
is the ratio g^|^o^^+p-j ■ The solid lines are drawn with central values of various parameters; the 
bands between the dashed lines denote measurement within la. The shaded dots originate from 
uncertainties due to variations of parameters, such as in Fig.(b) due to Xa, in Fig.(c) due to Xh, 
in Fig.(d) due to form factors, in Fig.(e) due to CKM elements, and in Fig.(a,f) due to various 
parameters. 
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FIG. 4. Br[B~ ^iT~ p^) versus 7 at the scale /i = m;, with the QCDF approach, including 
the weak annihilation contributions. The solid lines are drawn with the default values of various 
parameters; the bands between the dashed lines denote averaged measurement within la. The 
shaded dots demonstrate the uncertainties due to variations of various parameters in Fig. (a), Xa 
in Fig.(b), Xh in Fig.(c), form factors in Fig.(d), and CKM matrix elements in Fig.(e). 
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FIG. 5. Br{B~^7T~uj) versus 7 at the scale ^ = nif, with the QCDF approach, including the 
weak annihilation contributions. The meaning of the solid lines, the bands, and shaded dots is the 
same as in Fig.0. 
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FIG. 6. The CP-averaged branching ratios of B^PV versus 7 at the scale jj, = rrih with the 

QCDF approach, including the weak annihilation contributions. The meanings of the solid lines 

and bands are the same as in FigJ^. The shaded dots denote the uncertainties due to variations of 

various parameters. 
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FIG. 7. CP-violating asymmetries Acp{B^^K^(f>) in Fig.(a) and Acp{B^^Tr^u>) in Fig.(b) 
versus 7 at the scale /x = nifj with the QCDF approach, including the effects of weak annihilations. 
The meanings of the solid lines, bands, and shaded dots are the same as in Fig.0. 
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TABLES 

TABLE I. Wilson coefficients in NDR sclieme. Tlie input parameters in numerical calcula- 



tions arc fixed ; 


-IS Q.s{mz) 


= 0.1185, Qem 


{mw) = 1/128, 


mw = 80.42GeV, mz = 


: 91.188GeV, 


mt = 168.2GeV, rrib = 4.6GeV. 












M — 




/X = mb 




= 2m6 




1\TT 


T 


NLO 


LO 


NLO 


LO 


/-I 


i.iOD 


i.i (0 


1.080 


1.115 


1.044 


1.073 


/-I 


-U.Zoo 


-U.ooo 


-0.181 


-0.245 


-0.105 


-0.165 






U.UZU 


0.014 


0.012 


0.009 


0.008 


r~< 
G4 


-U.UoU 


-U.U4o 


-0.035 


-0.033 


-0.024 


-0.023 


r~< 


U.UiU 


U.UiU 


u.uuy 


U.UUo 


0.007 


0.006 


Ce 


-0.063 


-0.061 


-0.041 


-0.038 


-0.026 


-0.024 


Ct/ OLem 


-0.020 


-0.103 


-0.004 


-0.096 


n ni Q 
u.uiy 


-0.081 


C%/ aem 


0.082 


0.024 


0.052 


0.015 


0.033 


0.010 


Cg/aem 


-1.339 


-0.089 


-1.263 


-0.086 


-1.201 


-0.074 


Cw/Oiem 


0.369 


-0.022 


0.253 


-0.016 


0.168 


-0.012 


C-j^ 




-0.341 




-0.304 




-0.272 






-0.160 




-0.145 




-0.132 
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TABLE II. Numerical values of coefficients aij with the default values of various parameters 
for case I (the recoiled mesons Mi are vector mesons, and the emitted mesons M2 are pseudoscalar 
mesons) . 







iLt = mb/2 


/X = mb 


/X = 2mb 


ai,i 




1.077+0.033i 


1.055+0.0181 


1.038+0.0101 


a2,i 




-0.023-O.llOi 


0.018-0.0821 


0.057-0.0651 


03,7(10- 


-4\ 


91.649+44.722i 


70.951+24.2591 


49.938+14.0361 




-4^ 


-331.4-188.31i 


-299.26-152.671 


-270.97-129.361 


«2,/(io- 


-4\ 


-406.48-60.422i 


-355.33-57.1621 


-315.76-53.0651 


05,7(10- 


-4\ 


-99.038-56.834i 


-65.932-27.9341 


-39.693-14.6461 


«6,/(10- 


-4\ 


-558.53-169.85i 


-416.46-143.941 


-327.55-124.531 


«g,/(io- 


-4\ 


-600.85-40.033i 


-448.06-46.9881 


-352.8- 47.081 


«6,7^x(10- 


-4\ 


-484.0-147. 19i 


-483.15-166.991 


-474.92-180.561 


O6,7^x(10" 


-4\ 


-520.68-34.691i 


-519.81-54.5121 


-511.53-68.2631 


a7,/(10- 


-4\ 


0.539-F0.175i 


1.101+0.0861 


2.419+0.0481 


«I,7(10- 


-4\ 


7.383-0.3301 


4.759-0.8031 


2.934-1.1651 


ai,/(10- 


-4\ 


7.344-0.2091 


4.632-0.4151 


2.741-0.5741 


<7^x(10- 


-4\ 


6.398-0.2861 


5.521-0.9321 


4.254-1.6891 


«8,/''x(10" 


-4\ 


6.364-0.1811 


5.374-0.4811 


3.974-0.8321 


09,7(10- 


-4\ 


-94.827+0.1601 


-91.903+0.0921 


-89.29+0.0571 


aio,/(10- 


-4\ 


-2.426-h0.318i 


-9.989-0.3091 


-15.972-0.8121 


«fo,7(10- 


-4-) 


-2.497+0.4381 


-10.214+0.0731 


-16.314-0.2301 
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TABLE III. Summary of experimental data for the branching ratios (in units of 10~^) for 
several charmless hadronic B^PV . Inequality denotes 90% C.L. upper limits. The last column is 
for the averages of the uncorrelated measurements 



decay decay 



BaBar 



Belle 



CLEO 



Averaged 



B-^TT-p^ 



B -^TT UO 



B -^K oj 



B'^rjK*' 

B-^n-K*' 



B-^rjK* 



B-^R-p^ 



B-^K-uj 



B-^K- 



B^TT+K*- 



28.9±5.4±4.3 



< 39 



6.6tli±0.7 



:< 13) 



19.8^|^±1.7 



15.5±3.4±1.8 |T| 

22.11^^2^ ±3.3 
(< 33.9) 
<29 i 



(< 4) 
7.7lJ:^±0.8 



20.2i|^±3.3 



(< 35.7) 



11.2ti^±l. 



(< 28.^ 
< 9.4 1 



2i.2l:|:^±2.o 



1Q /I +4.2+2.1+3.5 
-'■^■^-3.9-2.1-6.8 



< 13.5 



< 10.5 



11.2i|^±1.4 



26.0±8.3±3.5 



27.6t^j±4.2 1T| 25.9±4.5 



10.4l|^±2.1 ^ 10.7±3.2 



ii.3ty±i.4 m 

i6.ot^:^±2.8 m 

(<32) 

io.oiy±i.4 m 

(< 21) 



5.4+|^±0.7 



(< 12.3) 

i3.8tti±i.6 m 



7.6^3:0^1-6 



(< 16) 



8.4t|0±l. 



(< 17) 
3.2t?;^±0.8 

(< 7.9) 
5.5t^:^±0.6 



8.0±1.7 
15.9±4.4 

7.5±2.8 

7.5±1.8 

18.0±3.2 
14.3±2.4 

24.5±7.1 

8.4±4.1 

1.8±1.1 

7.9±1.2 
23.8±6.1 
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TABLE IV. CP-averaged branching ratios (in units of 10~^) of decays B^PV for 6— xi tran- 
sitions with central values of various parameters. The results in columns 2~4 are calculated with 
A = 0.819, A = 0.2237, p = 0.218, and f) = 0.316, while the results in columes 5~7 are computed 
with A = 0.83, A = 0.222, p = 0.05, and f) = 0.381. 



decay 


NF 


QCDF 


NF 


QCDF 




modes 


Br 


Brf Brf+'' 


Br 


Brf Brf+^ 


Exp. 



=,0 



B°^7]uj 0.253 0.113 0.129 0.206 0.092 0.100 < 12 pT| 



B^^rj'cj) 0.0001 0.0003 0.0002 0.0001 0.0004 0.0003 < 31 |I| 



B-^TT-p^ 7.758 6.464 6.498 8.125 6.949 6.954 see Tab.ttll 



'^^K^K*^ 0.023 0.034 0.040 0.034 0.050 0.061 — 

0.132 0.198 0.210 0.191 0.270 0.293 — 

B^^K^K*^ — — 0.019 — — 0.019 — 

S°^7r-p+ 9.231 9.685 10.13 9.31 9.836 10.28 see Tab.|l] 



B^7r+p~ 21.59 22.72 23.36 19.8 20.70 21.34 see Tab. Ill 



S°^7rV 0.452 0.132 0.139 0.434 0.124 0.127 < 5.5 [|1 

B^^TT^iu 0.022 0.035 0.022 0.016 0.027 0.027 < 3 § 

S°^7rV 0.0003 0.0006 — 0.0004 0.0008 — < 5 El 



B^rip" 0.002 0.010 0.037 0.004 0.014 0.042 < 10 

B^^rj'p^ 0.009 0.019 0.040 0.007 0.016 0.030 < 12 |l| 



B^^r]'uj 0.161 0.061 0.068 0.162 0.063 0.072 < 60 |T| 



S°^?70 0.0001 0.0004 0.0005 0.0002 0.0006 0.0007 < 9 El 



B-^TT-uj 7.988 7.08 6.977 7.292 6.349 6.260 see Tab. Ill 



B-^Tr-(j) 0.0006 0.0012 — 0.0009 0.0017 — < 5 H 



B-^TT^p- 14.64 13.55 13.54 13.05 11.91 11.97 < 43 p| 



B-^'qp- 6.627 5.879 5.798 6.197 5.489 5.415 < 15 [|T3| 



B-^r]'p- 4.954 4.424 4.366 4.935 4.494 4.438 < 33 [|T3| 

B-^K^K*- 0.025 0.036 0.045 0.036 0.054 0.061 — 
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*o 



0.141 



0.20 



0.211 



0.204 



0.274 



0.311 



< 5.3 



TABLE V. CP-averaged branching ratios (in units of 10~^) of decays B^PV for 6— >s transi- 
tions with central values of various parameters. The results in columns 2~4 are calculated with 
A = 0.819, A = 0.2237, p = 0.218, and f) = 0.316, while the results in columns 5~7 are computed 
with A = 0.83, A = 0.222, p = 0.05, and f) = 0.381. 



decay 


NF 


QCDF 


NF 


QCDF 




modes 


Br 


Brf Brf^" 


Br 


Brf Br^^" 


Exp. 





^K-p+ 


1.485 


1.848 


2.008 


1.081 


1.321 


1.501 






0.918 


1.184 


1.256 


1.038 


1.239 


1.297 






0.034 


0.083 


0.007 


0.026 


0.076 


0.012 




.kV 


3.663 


5.945 


6.703 


3.589 


5.833 


6.569 


B - 




1.838 


2.411 


2.743 


3.281 


4.077 


4.36 


b'- 




0.533 


0.744 


0.896 


0.459 


0.714 


0.875 


b'- 




2.072 


2.681 


2.972 


2.15 


2.67 


2.927 


b'- 


-*r] K 


0.759 


1.717 


1.891 


0.689 


1.662 


1.84 


B-- 


-^vr K 


2.583 


3.497 


3.814 


2.531 


3.433 


3.731 


B-- 




1.852 


2.317 


2.489 


3.067 


3.543 


3.667 


B-- 




1.777 


2.247 


2.591 


2.564 


3.022 


3.306 


B-- 




1.446 


2.664 


2.83 


1.046 


2.091 


2.277 


B-- 


-^K p 


0.403 


0.598 


0.789 


0.395 


0.585 


0.777 


B-- 


-^K-p^ 


0.453 


0.426 


0.528 


0.609 


0.503 


0.631 


B-- 




0.583 


0.530 


0.435 


0.580 


0.565 


0.495 


B-- 


^R-6 


3.911 


6.346 


7.179 


3.831 


6.227 


7.02 



see Tab. ID 



< 39 |I| 
see Tab. ID 
see Tab. ID 
see Tab. Ill 

< 3.6 @ 
see Tab. ID 



< 24 [|T3| 
see Tab. ID 



< 31 PI 
see Tab.lnl 

< 35 I 
<48 |1| 

see Tab. ID 
see Tab.|Tl 
see Tab. ID 
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TABLE VI. CP-violating asymmetry parameters a^/, and Og+g/ (in units of 10~^) for decays 
B^^PV with central values of various parameters within the QCDF framework. The results in 
columns 2~5 are calculated with A = 0.819, A = 0.2237, p = 0.218, and fj = 0.316, while the 
results in columns 6~9 are calculated with A = 0.83, A = 0.222, p = 0.05, and fj = 0.381. 



modes 


4 


f+a 

ay 


f 






f+a 

ay 


f 

<+e' 


f+a 

4+e' 




-11.11 


-6.90 


48.29 


51.66 


-14.01 


-8.88 


-39.15 


-34.80 




19.08 


80.95 


94.80 


32.56 


28.82 


80.02 


95.35 


10.70 




31.78 


16.77 


-17.31 


10.53 


28.11 


17.11 


-87.88 


-75.06 




-28.02 


-34.78 


89.23 


92.18 


-39.96 


-54.92 


88.09 


66.99 


rjuj 


44.82 


29.66 


74.53 


81.05 


64.64 


45.06 


16.02 


26.05 


r)'LO 


-20.57 


-16.88 


32.27 


26.95 


-23.60 


-18.81 


-56.70 


-62.16 




-5.56 


-8.92 


64.39 


66.25 


-6.27 


-10.20 


62.38 


64.40 




-41.77 


86.90 


79.28 


-48.78 


-54.10 


62.18 


73.15 


-61.69 




-0.99 


-0.97 


73.24 


73.40 


-1.19 


-1.17 


72.91 


73.11 


ttV 







2.29 









1.83 












2.29 


2.29 








1.83 


1.83 
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TABLE VII. CP-violating asymmetries Acp (%) for B^PV {b^d transitions) with central 
values of various parameters within the QCDF approach. The results in the 2nd and 3rd columns 
are calculated with A = 0.819, A = 0.2237, p = 0.218, and f] = 0.316, while the results in the 4th 
and 5th columns are computed with A = 0.83, A = 0.222, p = 0.05, and fj = 0.381. 



modes 


Aqp 


A J +a 


A-' 


4 f-i-n 
Aj-ra 


B^^K^K*^ 


22.38 


23.89 


19.60 


20.37 


B^^K^K*^ 


-11.34 


-11.97 


-9.14 


-9.36 


B^^K^K*^ 




17.41 




-25.49 


i-»n j_ 


17.09 


20.16 


-13.47 


-9.60 


B^^TT+p- 


25.64 


21.51 


-26.91 


-31.85 


B -^TT^p^ 


15.75 


20.10 


-27.78 


-22.36 


B ^TT^LU 


57.59 


68.32 


64.21 


57.29 


B -^■qp^ 


12.49 


15.95 


-23.52 


-24.58 


B ^rfp^ 


24.21 


21.21 


15.88 


-3.92 




64.74 


57.95 


49.80 


41.81 


B -^r]'u 


1.95 


1.83 


-42.40 


-41.88 


B ^7r"0 


1.09 




0.87 




B ^r/W^ 


1.09 


1.09 


0.87 


0.87 


B-^TT-p^ 


3.26 


12.47 


3.58 


13.76 


B~^'rr~LO 


-6.22 


-6.64 


-8.18 


-8.73 


B-^7r-(j) 












B'^TT^p- 


-2.98 


-9.09 


-4.0 


-12.14 


B-^rip- 


-1.75 


-2.05 


-2.21 


-2.59 


B-^T]'p- 


5.03 


4.94 


5.85 


5.74 


B-^K^gK*- 


-8.82 


8.94 


-6.94 


7.83 


B-^K-K*^ 


-25.28 


-33.23 


-21.78 


-26.58 
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TABLE VIII. CP-violating asymmetries Acp (%) for B^PV {b^s transitions) with central 
values of various parameters within the QCDF approach. The results in the 2nd and 3rd columns 
are calculated with A = 0.819, A = 0.2237, p = 0.218, and f] = 0.316, while the results in the 4th 
and 5th columns are computed with A = 0.83, A = 0.222, p = 0.05, and fj = 0.381. 



modes 






f 

•^CP 


f 1 „ 

•^CP 


B^^K-p+ 


1.62 


-34.31 


2.68 


-54.21 


B^^K^^ 


27.04 


25.73 


25.62 


24.01 


B ^K^u 


10.51 


33.46 


-0.46 


11.19 


B ^K^s<P 


34.23 


34.32 


33.95 


34.05 




20.57 


47.0 


14.36 


34.92 


B -^TT^K 


-9.49 


-9.51 


-11.69 


-11.50 


B —>riK 


5.05 


5.52 


5.99 


6.61 


B -^r]'K 


-5.51 


-5.56 


-6.72 


-6.75 


B ^TT K 


0.97 


1.23 


1.17 


1.49 


B -^n K 


1 o or* 

18.86 


35.57 


14.56 


28.51 


B-^r]K*- 


7.73 


29.28 


6.79 


27.09 


B-^rj'K*- 


-13.78 


-20.39 


-20.73 


-29.92 


B-^K^SP- 


0.31 


-0.34 


0.38 


-0.41 


B-^K-p^ 


2.88 


-80.25 


2.87 


-79.31 


B-^K-uj 


59.85 


-19.92 


66.26 


-20.63 


B-^K-(j) 


0.99 


1.17 


1.19 


1.41 
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